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A model for parton distributions in hadrons* 
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The non-perturbative parton distributions in hadrons are derived from simple physical arguments resulting in 
an analytical expression for the valence parton distributions. The sea partons arise mainly from pions in hadronic 
fluctuations. The model gives new insights and a good description of structure function data. 



Hard processes involving hadrons are calcu- 
lated by folding perturbative QCD matrix ele- 
ments with parton distributions describing the 
probability of finding a quark or a gluon in the 
hadron. 

Perturbative QCD evolution describes the de- 
pendence of the parton distributions on the hard 
scale Q of the interaction. However, their depen- 
dence on the momentum fraction x at the lower 
limit for applying perturbative QCD, Q w 0.5—2 
GeV, are fitted to data using parameterisations, 
e.g. of the form [ g g | 

fi{x, Q ) = N iX ai (1 - x) b ^ (1 + Cl yfc + d lX ) (1) 

The parameters in these functions have no di- 
rect physical meaning, making it difficult to inter- 
pret the results. To gain understanding of non- 
perturbative QCD we have developed a physical 
model [||] for the parton distributions at Qq. This 
is here briefly described together with our latest 
developments. 

The basic physical picture is that a probe with 
large resolution, compared to the hadron size, 
will see free quarks and gluons in quantum fluc- 
tuations of the hadron. The measuring time is 
short compared to the life-time of the fluctuation 
since the latter is determined by the confinement 
of quarks and gluons inside the hadron, as illus- 
trated in Fig. |l|. This makes it possible to describe 
the formation of the fluctuations independently of 
the measuring process. 
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Figure 1. Fluctuation of a hadron into a parton 
and a remnant. 
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Figure 2. Kinematics for scattering on a parton 
in a hadron. 



Our approach only intends to provide the four- 
momentum k of a single probed parton. All other 
information in the hadron wave function is ne- 
glected, treating the other partons collectively as 
a remnant with four- momentum r, see Fig. g. 

It is convenient to describe the process in the 
hadron rest frame where there is no preferred di- 
rection and hence spherical symmetry. The prob- 
ability distribution for finding one parton is taken 
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as a Gaussian, which expressed in momentum 
space for a parton with four-momentum k and 
mass rtij is 
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fi[k)dk = N(ai, m,i)e 



dk, (2) 



where <r = -f- rs is the inverse of the confinc- 

d h 77 

ment length scale or hadron diameter. 

The partonic structure is described using the 

light-cone momentum fraction x — — which the 

. ... p+ 
parton has in the initial hadron. Since x is invari- 
ant under boosts in the z direction, the same will 
be true for the calculated parton distributions. 

There are a number of constraints that must 
be fullfilled by the parton distributions. The nor- 
malisation for valence quarks is given by the sum 
rules 



fi(x)dx = Tli, 



(3) 



and for the gluons by the momentum sum rule 



y~] / xfi(x)dx = 1. 
. Jo 



(4) 
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Figure 3. Valence distributions at Qo 
in the proton and tt + . 
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There are also the kinematical constraints 



< j 2 < W 2 



and 



r 2 > 



(5) 



given by the final partons being on-shell or time- 
like and the remnant having to include the re- 
maining partons. These constraints also leads to 
< x < 1. 

The parton model requires that W is well above 
the resonance region and that the resolution of 
the probe is much larger than the size of the 
hadron, i.e. 



W > m p and Qo > cr; 



(6) 



The scale of the probe must also be large 
enough, Qq 3> Aqcd, for perturbative QCD to 
describe the evolution of the parton distributions 
from the starting scale Qq. 

In [ ||] we integrated Eq. [2] numerically to find 
the parton distributions since the kinematical 
constraints, Eqs. [| are quite complicated in gen- 
eral. The problem is much simpler if the trans- 
verse momenta and the masses of the partons are 



neglected. It is then possible to derive [ |) an an- 
alytical expression for the parton distributions, 



/;(. r) : V'(o u xpl-^erf 



where 
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(7) 



(8) 



The valence parton distributions for hadrons 
are here determined simply by the mass and size 
of the hadron! The resulting valence distribu- 
tions for the proton (a ~ 0.15) and the pion 
(cf w 1) are very reasonable as shown in Fig. ||. 
Note that the pion distributions are very similar 
to xf(x) = 2x(l — x) and that one third of the 
pion momentum is carried by gluons. 

The sea partons are described by hadronic fluc- 
tuations, e.g. for the proton \pir a ) + \nir + ) + 
where the probe measures a valence parton in one 
of the two hadrons. The momentum distribution 
of pions in a hadronic fluctuation is assumed to 
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Figure 4. Data from NMC [ [| and ZEUS [ 0] on the proton structure function F2(x,Q 2 ) compared to 
the model with valence partons only (dashed) and including sea partons (full). 



follow from the same model as for the valence 
partons, with the differences that the mass can 
not be neglected and the width ov ss 50 MeV 
is smaller, related to the longer range of pionic 
strong interactions. 

Using these valence and sea parton x distribut- 
ions at Qo = 0.85 GeV, next-to-leading or- 
der DGLAP evolution in the CTEQ program 
[ was applied to obtain the parton distribu- 
tions at larger Q. The proton structure function 
F2(x,Q 2 ) can then be calculated and compared 
with deep inelastic scattering data, as illustrated 
in Fig. |] and detailed in [ Q . The model does 
remarkably well, in view of its simplicity and few 
parameters (6). Of course, conventional parton 
density parameterisations give much better fits, 
probably mainly due to their many more param- 
eters (w 20). 

The main part of the proton structure is deter- 
mined by the valence distributions, but the sea 
gives an important contribution at small x, as 
can also be seen from the comparison with the 
measured F2. 

Including all pion fluctuations will give a flav- 
our asymmetric sea with d > u as also observed 
experimentally [ |), but the numerical details re- 
main to be investigated. 

The model predicts the valence parton distri- 



butions for all hadrons, but heavy quarks gives 
more complicated analytical expressions [ || . Nu- 
merical results on strange and charmed mesons 
are shown in [ Q . In addition, a study based on 
Monte Carlo has been made [ ^| to investigate 
intrinsic strange and charm quarks in the proton. 
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